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We discuss a possible new resonance in active-sterile neutrino oscillations arising in 
theories with large extra dimensions. Fluctuations in the brane effectively increase the 
path-length of active neutrinos relative to the path-length of sterile neutrinos through 



the extra-dimensional bulk. Well below the resonance, the standard oscillation for- 
mulas apply. Well above the resonance, active-sterile oscillations are suppressed. We 
show that a resonance energy in the range of 30 — 400 MeV allows an explanation of all 
lO ■ neutrino oscillation data, including LSND data, in a consistent four-neutrino model. 

o ■ 

A high resonance energy implies an enhanced signal in MiniBooNE. A low resonance 
energy implies a distorted energy spectrum in LSND, and an enhanced depletion 
from a stopped-pion source. The numerical value of the resonance energy may be 
related back to the geometric aspects of the brane world. Some astrophysical and 
cosmological consequences of the brane-bulk resonance are briefly sketched. 



1 Introduction 

Theories with large extra dimensions typically confine the Standard Model (SM) particles on 
a 3+1 brane embedded in an extra-dimensional bulk [UI2] see a ^ so 0- Gauge singlet parti- 
cles may travel on or off the brane. These include the graviton, and any singlet (= "sterile") 
neutrinos jlj. Virtual gravitons, too, penetrate the bulk, and so lead via Gauss' Law to an 
apparent weak gravity on our brane, when in fact the strength of gravity may unify with 
the SM forces at scales as low as a few TeV. Here we focus on the "other" possible particle 
in the bulk, the sterile neutrino [5]. 

The only evidence to date for the existence of the sterile neutrino comes from the 
incompatibility of all reported neutrino oscillation results with the three-active neutrino 
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world. The solar and atmospheric disappearance data are corroborated, whereas the LSND 
appearance data |Jj are not yet corroborated. We will assume that the LSND data is correct, 
and use this as motivation to study the possible compatibility of all the data when the higher 
dimensional bulk is included. We find a new active-sterile resonance which relates bulk and 
brane travel times. When the new resonance energy falls between the LSND energies and 
the CDHS energies, then all the oscillation data become compatible. 

This article is organized as follows: In section 2 we discuss a metric with small-scale 
fluctuations, which allows for bulk shortcuts. In section 3 the effect of bulk shortcuts on 
active-sterile neutrino oscillations is illustrated in a simple model with one sterile and one 
active neutrino. Section 4 discusses the LSND result in the context of constraints from 
other neutrino experiments in a realistic 3+1 neutrino model. Some further astrophysical 
and phenomenological constraints are discussed in section 5, and conclusions are drawn in 
section 6. 



It has been shown that branes embedded in higher dimensional spacetime are curved extrin- 
sically by self-gravity effects in the presence of matter [7j. For example, while the on-brane 
distance between atomic constituents is fixed by electromagnet ism and the Pauli principle, 
the attractive force of gravity between constituents can shorten the embedding distance of 
the brane in the bulk, leading naturally to a scenario where the brane is deformed (possesses 
fluctuations or "buckles") on a microscopic scale. Alternative causes of brane bending in- 
clude thermal and quantum fluctuations. This picture leads to a framework in which the 
on-brane geodesic felt by an active neutrino is longer than the bulk geodesic felt by a sterile 
neutrino. Such apparent superluminal behavior for gauge-singlet quanta has been noted 
before, for the graviton [El El- 
ba the spirit of 9] it is straightforward to construct a 1+1 dimensional toy model with a 
metric which indeed exhibits the anticipated behavior. First, write down a 1+2 dimensional 
embedding spacetime, with Minkowski metric 



In this embedding spacetime, assume the brane exhibits periodic (for simplicity) oscillations 
in space 



2 A metric for bulk shortcuts 



ds 2 = dt 2 — dx 2 — dy 2 . 




y = A sin kx ; 



(2) 
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Figure 1: Schematic representation of a periodically curved brane in Minkowski spacetime. 
A coordinate transformation leads to an equivalent description as a non-diagonal metric 
with a flat brane, as described in the text. 

see Fig. ^ for a schematic representation. Here k is the wave number of the fluctuation (in 
the x-direction), and A is the amplitude of the fluctuation (in the ^-direction). A brane 
with tension is dynamical, and fluctuations change in time. The toy-model fluctuations 
should be thought of as some rms average over many time-slices. 

The bulk-geodesic for the sterile neutrino is simply given by y = 0, which leads to a 
travel distance of 

Dg = x . (3) 
The geodesic for the active state on the brane is slightly more complicated: 

D b = f ^dx 2 + dy 2 = f X Vl + A 2 k 2 cos 2 kx dx . (4) 

J brane J 

We use subscripts b and g to denote the brane and bulk spaces, respectively. 

In terms of the coordinate x, the parameter describing the shortcut in the bulk is 

D b r vi + A 2 k 2 cos kcc dx 

While mathematically correct, this description of the geodesies as functions of x has a 
shortcoming, in that x is not a coordinate easily identified in an experiment on the brane. 
It is useful to consider a more physical set of brane coordinates. They will lead to essentially 
the same parameter e. 

Consider the space-coordinate transformation 

u = y — A sin kx (6) 



and 



Under this transformation, the line element in (JTJ transforms into 

o 7 2 t 2 72 *2 A k cos kx(z^) 
as — at — dz — au = du dz. 



+ A 2 k 2 cos 2 /co^z) 

Note that in (u, z) coordinates, u = defines the location of the brane, and z labels the 
physical distance along the brane. Consequently, a photon moving along the brane (u = 0) 
satisfies the equation 

ds 2 = dt 2 - dz 2 , (9) 
and thus travels in time tf the distance 

z b = t f . (10) 

On the other hand, in (x, y) coordinates, the brane is described by the periodic sine 
function, while the geodesic in the bulk follows a straight line along mean y, given by 

y 9 = o, D g = t. (li) 

Using eqs. (jUJ) and (J7J), the bulk geodesic equation (fTTj) can be transformed into the (w, z) 

system, 

u g = -A sin(H), (12) 
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C vi + ^wh dt = *1+±L £ ktf , J^- 2 , (is) 



where S(p, q) denotes the elliptic integral of the second kind. 

The bulk geodesic intersects the brane at u g = 0, which according to eq. (fT2"j) occurs at 
the discrete times 

tint = -y, (14) 

where n is an integer. However, if the size of the brane's fluctuations is small on the scale of 
an experimental detector, then it is not required that the two geodesies be in intersection, 
and t int has no special significance. 

From a comparison of the integrand in (|13|) to the result of (|l(Jj). one readily infers 
that z g > zj), which means that in a common time interval the bulk test particle seemingly 
travels farther in the physical z-coordinate than the brane particle. In other words, the 
specific metric (0) allows apparent superluminal propagation. The shortcut in the bulk can 
be parametrized by 



We note that this e and the one defined in eq. (j5J) are formally the same when the space- 
coordinate x is replaced by the physical time-coordinate t. 

The parameter e depends very weakly on tf when many fluctuations are traversed, i.e. 
when tf ^> 2ir/k. In fact, we are free to choose tf according to eq. ([14)1 . With this choice, 
the shortcut parameter depends only on the geometry Ak of the brane fluctuation, according 
to 

tt/2 



A 2 k 2 



(16) 



(1 + A 2 k 2 ) S c (J, y i+^fe 2 
In the latter equation, we have used the relation £(nir,q) = 2n£ c (|,g); the expression 



S c (f is called the complete elliptic integral. In the model developed here, an inference 
of e offers a direct measurement of the brane-fluctuation shape-parameter Ak. 

The dimensionless ratio Ak is essentially the aspect ratio (height to width) of the 
fluctuation. For a brane with high (low) tension, we expect a low (higher) value of Ak. 
Our assumption is that the brane tension is large, so that there is no curvature on large 
scales. Accordingly, we expect a small value for Ak. To first non-vanishing order in Ak, the 
parameter e is 

-(f) 2 - 

This approximation is valid until Ak becomes of order unity, after which z g 3> Zb and e 
itself approaches unity. 

An alternative to the periodic metric described here arises in spacetimes in which the 
speed of light along flat 4D sections varies with the extra dimension coordinate u due to 
different warp factors for the space and the time coordinates ("asymmetrically warped" 
spacetimes) [3 QDJ. Such scenarios are realized, e.g. by a black- hole background in the 
bulk, and may provide interesting consequences for the adjustment of the cosmological 
constant [ID] . A specific example is the two-brane scenario discussed in [DJ, in which the 
4+1 dimensional metric is given by 

ds 2 = dt 2 - [e~ 2ku a 2 (t) dh 2 + du 2 }. (18) 

Here h denotes an Euclidian three- vector. In this model a sterile neutrino could scatter out 
of our brane at u\ on a geodesic perpendicular to h, reach a second brane at u-x > Ui, and 
finally scatter off some fields confined on this second brane to return back to our brane. It 
has been shown in 19] that again the on-brane distance traveled in a given time interval via 
a path on the hidden-sector brane and in the bulk can be larger than the distance for pure 
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on-brane travel in the same time interval. We anticipate that the qualitative features of 
this and similar scenarios can be approximated by the simple toy model discussed above. 



3 The two-state sterile-active oscillation probability 

Let us illustrate the brane-bulk resonance for a simple system of one sterile neutrino v s and 
one active neutrino v a . The mass eigenstates are m,2 and mi, respectively, in the sense that 
for small mixing u s is mostly v-i and v a is mostly v\. 

Ignoring the bulk for the moment, the evolution equation in flavor space reads 

(19) 

and the Hamiltonian in the flavor basis is 

„ 1 5m 2 ( - cos 29 sin2# \ / s 

H F = E + -Tr + —-\ 20 
2 4E y S i n 26 cos2# ) 

where Tr = m\ + m\/2E, 5m 2 = m\ — m\ and 9 is the mixing-angle in the unitary matrix 
relating flavor and mass bases: 



= U* aj \uj) , or Uaj = (v a \Vj) , (21) 



with 



, cos 9 sin 9 , 
U=[ | ■ (22) 

\ —sm9 cos 9 

We will call 9 and Sm 2 the standard mixing angle, and standard mass-squared difference, 
respectively. They are the analogs of vacuum values in MSW [llj physics. 

Now we let the sterile neutrino propagate in the bulk as well as on the brane. If the 
brane were rigid and flat in its embedding, then the sterile geodesic is just the same as 
the active geodesic on the brane. However, if the brane is curved in its embedding, as 
discussed above, then the sterile neutrino may have a different trajectory, with a shorter 
geodesic than that of the active neutrino constrained to the brane. We will formulate this 
as an effective potential contributing to the sterile-sterile term of the Hamiltonian in flavor 
space. Note that this is analogous to the Wolfenstein potential for the active-active term 
due to forward elastic scattering in matter, albeit with three important differences. The 
first is that the effective potential here is the same for neutrino and antineutrino, because 
gravitationally-determined geodesies do not distinguish between particle and antiparticle. 
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The second difference is a more pronounced energy dependence here, with the effective mass- 
squared difference varying as E 2 , not as E. The third difference is that there is no time 
or space dependence in the Hamiltonian. We note that characteristics of this brane-bulk 
luminal/superluminal system resemble certain scenarios with Lorentz invariance violation 



Due to the shortcut in the bulk (see the schematic in Fig. (jTfl). the sterile state will 
appear to cover more distance on the brane than the active neutrino does in the same time, 
or equivalently, the same distance but in a shorter time. The ratio of apparent times at 
common distance or apparent distances at common time for the sterile and active neutrinos 
is dt/t ~ 5z/z ~ e. For the toy model introduced in section |21 the parameter e is given by 
eq. (|TB|) . or by eq. (fTTj) to lowest non- vanishing order in Ak. 

Adding the new contribution to the sterile-sterile element of the Hp, action/time= 
£f , and then removing the irrelevant energy and trace terms, we arrive at the effective 
Hamiltonian 



The bulk term may beat against the brane term to give resonant mixing, i.e., for some 
energy E Tes even a small standard angle can become large or even maximal in the brane- 
bulk model. The resonance condition is that the two diagonal elements in Hp be equal, 
which implies 



Since the value of e is unknown, the resonance energy could have almost any value, a 
priori. However, if e <C 1, as we assume, then we have the result 5m 2 -C E 2 es . Still, there 
is much parameter space available for resonance. Our aim is to accommodate the LSND 
result in a four-neutrino framework, and so we will restrict the resonance energy with this 
in mind. It is worth noting that according to (j2U), a determination of E ies fixes e, if 5m 2 
and cos 26 can be independently determined. One way to independently determine 5m 2 
and cos 26 is to observe the active-sterile oscillation parameters far below resonance, where 
the oscillations are described by the standard formulas. Note that knowledge of e, when 
available, yields the shape-parameter Ak of the brane fluctuation, according to eq. (fTTj) . 

The value of E ves naturally divides the energy domain into three regions. Below the 
resonance, oscillation parameters reduce to their standard values and give the familiar 
oscillation results. At resonance, the mixing angle attains a maximum (but the effect on 
the oscillation probability can be reduced by a compensating factor in the 5m 2 term). Above 



H2i unl- 




ess) 




(24) 
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resonance, the oscillations are suppressed. Our strategy to accommodate the LSND data in 
a four-neutrino framework will be to set the resonant energy well below the CDHS data to 
suppress oscillations for this experiment, but at or above the LSND energies, so as to not 
suppress (or even, to enhance) the LSND signal. 

To find the new eigenvalue difference 5H and the new mixing angle 9 effected by the 
bulk, one diagonalizes the 2x2 system. In terms of the new 5H and 9 one obtains the 
usual expression for the flavor-oscillation probability 



P as = sin 2 29 sin 2 (SHD/2) 

with new values given in terms of standard values by 

2 n sin 2 20 
sin 2 29 = 



(25) 



5H 



sin 2 29 + cos 2 29 
5m 2 



1 - 



E 



2 ■ 



2E \ 



sin 2 29 + cos 2 29 



E 



21 2 



(26) 



(27) 



The width of the resonance is easily derived from the classical amplitude in eq. pB*]) . 1 Fig. [2] 
shows sin 2 29 for different values of sin 2 29 as a function of energy. 

For E 3> E res , the sterile state decouples from the active state, as sin 2 29 — > 0. Although 
the example presented in this section contains a single sterile state and a single active state, 
the decoupling of the sterile state (s) from the active state (s) is a general feature. 



4 Accommodating the LSND result 



As the sterile neutrino mass is not protected by the gauge symmetry of the Standard Model, 
it is natural to assume it to be larger than the masses of the active neutrinos. We thus 
focus here on a 3 + 1 neutrino spectrum [T%] . i.e. three active neutrinos are separated by 



the LSND mass-squared gap 5^ 2 snd from the dominantly sterile state z/ 4 ~ v s . 
1 A short calculation gives the Full Width in energy at a fraction / of Maximum (FWfM) as 

1/2 r , -. 1/2 



A£(FWfM) 



1 + tan 29 , 



1 - tan 26* , 



(28) 



which, for small 9 reduces to 29 y ~f- Thus, the resonance is very narrow for a small standard angle. For 
example, Full Width at Half Max is AiJ(FWHM) = 2 9 E rcs for small angle. For a larger standard angle, 
the resonance becomes less dramatic. 
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Figure 2: Oscillation amplitude sin 2 26 as a function of the neutrino energy E u , for a 
resonance energy of E ieB = 40 MeV. The different curves correspond to different values for 
the standard angle, sin 2 26 = 0.2,0.1,0.01,0.001 (from above). 





a(3 


E v 


D 


sin 2 29 a/ 3 


LSND 


fie 


20-52.8 MeV 


30 m 


> 0.003 


KARMEN 


fie 


20-52.8 MeV 


17.7 m 


< 0.002 


MiniBooNE 


fie 


0.1-1.0 GeV 


540 m 


~ 0.0006 


BUGEY 


ee 1 


1-6 MeV 


25 m 


< 0.15 


CDHS 


//// 


> 1 GeV 


755 m 


< 0.1 



Table 1: Flavor channels, beam energies, oscillation distances (or differences of far and 
near detectors), and limiting oscillation amplitudes in the large 5m 2 limit, for the relevant 
experiments [HI El E3 E] ■ For BUGEY, the 25 m measurement has been chosen. 
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In the present model, (^lsnd differs from the standard formalism by the square-root 
factor in (|27j) . If the new resonance occurs in the energy range of LSND/KARMEN, this 
factor might allow a fit with a larger 5m 2 . We do not pursue this subtlety in this work (see, 
however, the discussion in section I5~2*j) . 

For short distances, the mass-squared differences in the 3 + 1 spectrum can be taken 
as Sm 2 = 5m 2 4 = Sm 2 ^ = Sm 2 A , with all other differences set to zero. There results in all 
short- distance oscillation probabilities a universal factor of sin 2 (8H D/2), with 5H given 
previously in eq. (J27|) . As a result of this universality, the relevant oscillation amplitudes 
can be completely described using a two-neutrino formulation. First, we define the linear 
combination of active flavors which couples to the heavy lz/4) mass eigenstate as |z/ a ). Below 
the resonant energy, we write 

lz/4) = cos 9\v s ) + sin #|z/ a ) . (29) 

Above the resonant energy, the mass eigenstates have to be relabeled (4 <-» 1); i.e. the 
isolated state contains little \v a ) but much \v s ). Put another way, the oscillations considered 
still occur above the resonance over the large active-sterile mass gap, but effectively with the 
interchange cos 9 <-> sm9 in eq. (J29)) . We note that sin 2 29 is unaffected by this interchange. 

To produce z/ M -z/ e oscillations, it is necessary that this state |z/ a ) contains |z/ M ) and |z/ e ). 
It may also contain |z/ T ). For simplicity, we will take this state to be a mixture of just |z/ M ) 
and \v e ): 

\u a ) = cos#*|^) +sin0*|i/ e ) . (30) 

Thus, we have 

[/ e 4 = sin 9 sin 9, , (31) 
UpA = sin 9 cos 9, , (32) 

below the resonant energy, and cos 9 <-> sin 9 in eqs. (JBTj) and (J3*2~)) above the resonant energy. 

The oscillation amplitude relevant for the LSND [U], KARMEN [Hj and MiniBooNE 
[To] appearance experiments is given by 

sin 2 2^ LSND = -4 ]T U eA U ej % 4 = 4 U 2 eA U 2 4 = sin 2 29, sin 4 9 . (33) 

j<4 

Similarly, the oscillation amplitudes for the v e and disappearance experiments BUGEY 
[To] and CDHS [T7] are given by 

sin 2 29 V = 4 U 2 4 ( 1 - U 2 4 ) = 4 sin 2 9 sin 2 9,(1- sin 2 9 sin 2 9* ) , (34) 
sin 2 29^ = 4t/ 2 4 (l - U 2 ^) = 4 sin 2 ^cos 2 ^ (1 - sin 2 £ cos 2 ft,), (35) 
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respectively. The far-right expressions in these equations hold below the resonance; above 
the resonance one must interchange cos 9 «-> sin 6, as explained above. 

Eq. (J53jl may be inverted to give £/ 2 4 = sin 2 9 e ^, and eq. to give f/ 2 4 = sin 2 6^. 
Were these probabilities not energy-dependent, as in the standard vacuum case, one could 
substitute these results into eq. to get 

sin 2 2#lsnd = 4 sin 2 6^sin 2 6^ ~ ^ sin 2 2# e ^sin 2 29^ , (36) 

with the latter expression holding in the small angle approximation. One recovers the 
well-known result in the standard case, that the LSND amplitude is doubly suppressed by 
stringent bounds on the BUGEY and CDHS amplitudes. This fact excludes the standard 
3+1 neutrino models from describing the results of all short-baseline neutrino experiments 
[TSj . However, eq. ()36j) is not valid in general in our bulk shortcut scenario. As we shall 
demonstrate, the energy- dependence imparted to the mixing-angles and to 5H by the bulk 
shortcut allows consistency. 

It will also be useful to list the amplitudes for u^-Ug and v e -v s oscillations, for these 
may affect atmospheric and solar oscillations. As shown in eqs. ()26I27|) . the active-sterile 
oscillations are governed by vacuum values below the resonant energy, are maximal at the 
resonant energy, and are suppressed above the resonant energy. The same is true therefore 
for v^-v,, and v e -v s oscillations: 

sin 2 29^ = cos 2 0* sin 2 26 , (37) 

and 

sin 2 29 es = sin 2 6* sin 2 29 , (38) 

with sin 2 26 given in eq. (|26|) . 

For the analysis of neutrino oscillations in the bulk-shortcut scenario, the energy of the 
neutrino beam is of crucial importance. In Table Q the relevant experiments are shown 
together with the flavor amplitude to which they are sensitive, the neutrino beam energy, 
and the bound (or for LSND, the favored region) for the amplitude. Our task is to compare 
these experimental bounds (and for LSND, the positive signal) with the energy-dependent 
oscillation probabilities listed above. Note that all experimental data must be accommo- 
dated with four parameters: the standard mass-squared difference 5m 2 , the standard mixing 
angles 9 (describing v a -v s mixing) and (parametrizing the flavor composition of u a ), and 
the shortcut parameter e (or equivalently, the resonant energy En, as given in eq. (|24j) ) . 
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4.1 BUGEY 



BUGEY detects reactor neutrinos of energies an order of magnitude below the energies of 
LSND. We assume that the BUGEY energies are also far below the resonance energy, in 
which case the bulk shortcut effects decouple and the oscillation amplitude is given by the 
standard limit of eq. (JsHJ), 

sin 2 26ej ~ 4 U 2 eA ~ 4 sin 2 9 sin 2 6* . (39) 

Here the smallness of sin 2 29 ef / has been assumed, in concordance with the BUGEY limit 
sin 2 29 ef ^ < 0.15. We note that the BUGEY amplitude maybe suppressed by a small sin 2 
or small sin 2 9, or both. 

In the following we discuss two parameterizations compatible with LSND and CDHS 
which sufficiently suppress this BUGEY amplitude, namely sin 2 29 = 0.9, sin 2 9* = 0.01, and 
sin 2 29 = 0.45, sin 2 9* = 0.1. The first parameterization has a large active-sterile mixing, 
9 = 36°, but still a small sin 2 29 e ^ = 0.014. The second parameterization yields a moderate 
9 = 21° and small sin 2 29 efj t = 0.052. Since the BUGEY amplitude is sufficiently suppressed, 
there is no bound from BUGEY data on the value of 8m 2 . 

Future reactor experiments are proposed to search for nonzero U e ^. In these experiments, 
mixing of the sterile state with v e would have a measurable effect, mimicking a non-zero 
U e s. In contrast to the effects of a non-zero U e ^ (equivalently, a non-zero #13) and as a result 
of the large oscillation phase oc Sm^^j), the effect will be seen in both the near and the far 
detectors, though, and can be as large as sin 2 29 ri = 0.05. 

4.2 CDHS 

The accelerator oscillation experiment CDHS operated with neutrino energies above a GeV. 
At energies E 3> E res the active-sterile mixing is suppressed, and one can approximate for 
small sin0* (from eqs. (pjj) and (I3*B1 ) 

sin 2 29^ ~ cos 2 9* sin 2 29 ~ cos 2 9* tan 2 29 ( — — J . (40) 

This implies that neutrino oscillations in the CDHS experiment are suppressed by a factor 
between 10 6 and 40 for a resonance energy in the range of 30-400 MeV, making the oscillation 
amplitude sin 2 29^ unobservable above ~ GeV, even if the sin 2 29^ were maximal below 
the resonance. 
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4.3 LSND and KARMEN 

We have suppressed the BUGEY oscillation amplitude with the choice of a small U 2 A below 
resonance. We have suppressed the CDHS oscillation amplitude with the choice of a resonant 
energy below 400 MeV. This leaves two possibilities for the effect of the resonance on the 
LSND/KARMEN energy range, 20-53 MeV. The resonance may occur within this range, in 
which case the effect is observable. Alternatively, the resonance may occur above this range, 
in which case there is little change from the standard prediction and fit. We show below 
that either possibility can be realized with a resonance consistent with all other oscillation 
data. We also show below that the two different choices have very different consequences 
for the on-going MiniBooNE experiment. 

From eq. (J3*3j) and the following discussion, we have for the LSND oscillation amplitude 

sin 2 2fl LSND = X - sin 2 29* ( 1 =f cos 26) 2 . (41) 

The sign of the cos 29 term corresponds to energies below and above E res , respectively, 
accounting for the fact that the states have to be relabeled when crossing the resonance so 
that the oscillations considered occur over the large active-sterile mass gap. This formula 
also applies for the KARMEN and MiniBooNE amplitudes. 

In Figs. Q and (£Q) the oscillation probability P^ e predicted for LSND in the bulk- 
shortcut scenario is compared to the prediction of the standard oscillation case (dashed). 
Also shown is the expectation for KARMEN. Two very different sets of parameters have 
been chosen: one scenario has a resonant energy E res = 33 MeV in the LSND/KARMEN en- 
ergy range, and the other has a resonant energy E ves = 400 MeV far beyond the LSND/KARMEN 
energy range. 

To maintain consistency between the LSND and KARMEN data, it is necessary to 
exploit the differing distances of the two experimental configurations, £>lsnd = 30 m and 
-Dkarmen = 17.7 m according to Table 1. As in the standard approach, this is done as 
follows: According to eq. (|25|). the neutrino remains in its first oscillation until 5H D = 
2 n. For 5HD < 2n, the factor sin 2 (5H D/2) is well approximated by just (5H D/2) 2 , 
giving oscillation probabilities a quadratic dependence on distance. Since the baseline for 
KARMEN is about half that of LSND, choosing parameters such that 5H -Dkarmen ^ 1, 
with SH given in eq. (J2ZJ), suppresses KARMEN by a factor of four compared to LSND. 
This allows a slice of LSND parameter space to remain viable, in the face of the KARMEN 
null result. 

The requirement for LSND/KARMEN neutrinos to remain within their first oscillation 
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Figure 3: Bulk shortcut oscillation probabilities for LSND (light /green) and KARMEN 
(dark/red) as a function of the neutrino energy. Shown is a scenario with E Tes = 33 MeV; 
sin 2 0* = 0.01; sin 2 28 = 0.9; 5m 2 = 0.7 eV 2 . For comparison, a standard oscillation 
probability for LSND (5m 2 = 0.8 eV, sin 2 2fl LSND = 0.006) is displayed (dashed). The 
vertical lines indicate the energy window of LSND and KARMEN. 

is the standard one, 5m 2 ~ eV 2 . We adopt this value here. 

At this point, we may invert eq. (|24|) to determine the value of e: 



The freedom for E Tes in this model allows e to range over ~ 10~ 18 to 10~ 16 . According 
to eq. (|17|). this in turn implies a shape-parameter (height to width ratio) for the brane 
fluctuation of Ak ~ 10 -8 . These parameter values have to be eventually explained in a 
theory of brane dynamics. 

Both choices for E Tes exhibit a viable LSND energy spectrum, as seen in the figures. 
For the case where E vcs 3> -Elsnd, the LSND/KARMEN analysis and fit is the standard 
one. For the case where E res lies in the LSND range 20 MeV< E < 53 MeV, the energy 
dependence of the oscillation amplitude is modified considerably. For this latter case, we 
expect that the resonance should be evident in the LSND spectral data, and we encourage 
a re-analysis of the measured LSND energy-spectrum by the collaboration. 




(42) 
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Figure 4: Bulk shortcut oscillation probabilities for LSND (light/green) and KARMEN 
(dark/red) as a function of the neutrino energy. Shown is a scenario with E ies = 400 MeV, 
sin 2 6* = 0.1; sin 2 28 = 0.45; 5m? = 0.8 eV 2 . For comparison, a standard oscillation 
probability {5m 2 = 0.8 eV, sin 2 2# LS nd = 0.006) for LSND is displayed (dashed). The 
vertical lines indicate the energy window of LSND and KARMEN. 

4.4 Solar and Atmospheric Data 

The active-sterile amplitudes for and v e are given in eqs. (J37j) and (|3*5j). However, the 
same physics which suppresses active-sterile oscillations in BUGEY data and CDHS data, 
also suppresses active-sterile oscillations in solar and atmospheric data, respectively. For 
solar oscillations, the smallness of \U e ^\ 2 evades the experimental constraint. For atmo- 
spheric neutrinos, the measured energies are above the resonant energy, and so atmospheric 
oscillations into the sterile state are suppressed. The event sample below 500 MeV may 
contain enhanced v s production, but experimentally this would be difficult to confirm. 

In fact, even the 2+2 model of four-neutrinos can be resurrected with the brane-bulk 
resonance. In the same way that the suppression of sterile-active oscillations above the 
resonant energy neutralizes the CDHS constraint for the 3+1 model, so does it neutralize 
the atmospheric constraint for the 2+2 model. 

A detailed calculation is required to determine the nearly unitary 3x3 active-neutrino 
mixing-matrix that results when the sterile state decouples at high energy. We do not 
pursue this here. However, we expect that the freedom to partition the state \ u a ) among the 
\u e ), {v/j), and \u T ) is sufficient to yield acceptable phenomenology. For example, although 
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Figure 5: Bulk shortcut oscillation probabilities for MiniBooNE as a function of the neu- 
trino energy. Shown is a scenario with sin 2 9* = 0.1; sin 2 29 = 0.45; 5m 2 = 0.8 eV 2 . 
The resonance energy is varied, E res = 200, 300, 400 MeV, from left to right (light to 
dark). For comparison, the expectation for a standard oscillation solution (5m 2 = 0.8 eV, 
sin 2 2#lsnd = 0.006) for LSND is displayed (dashed). The vertical line indicates the energy 
threshold of MiniBooNE. 

we have taken (u T \u a ) = for simplicity, a «-> v T interchange symmetry, known to be 
consistent with all present data, can be incorporated here by changing \v^) in eq. (|3U|) to 



The requirement that E res < 400 MeV, well below the CDHS energy of ~ 1 GeV, leaves the 
resonance energy in the MiniBooNE range, 0.1 to 1 GeV, or even below. We predict that 
MiniBooNE should see no signal above ~ 700 MeV in this model. 

If E TCS falls in the MiniBooNE range above 100 MeV, then MiniBooNE should observe 
a strongly enhanced signal as evidence for the bulk-shortcut resonance. Near the resonance 
region, the exact expression (j4*Tj) applies for MiniBooNE. Results of this expression are 
shown in Fig. 03 for resonance energies of 200, 300 and 400 MeV. As can be seen, the 
strongly enhanced oscillation probability is unmistakable. 

On the other hand, if E res lies below the MiniBooNE threshold energy, then active-sterile 
mixing is strongly suppressed for MiniBooNE. At energies E ^> E res , one uses eq. (|2*Bj) in 
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4.5 MiniBooNE 
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Thus a null result is predicted for MiniBooNE in the case of a resonance energy (as in our 
33 MeV example) below the MiniBooNE threshold of 0(100) MeV. 

However, if E res is too low for an observable effect in MiniBooNE, a distortion in the 
LSND spectrum is expected (see Fig. (2J) • A strong disappearance signal also is predicted 
for an experiment using neutrinos from stopped pions, being proposed for the Spallation 
Neutron Source (SNS). This we discuss next. 

4.6 Muon-neutrino disappearance at the SNS 

While the sterile neutrino effectively decouples from the active sector at the CDHS energy 
and above, there is no suppression of the active-sterile mixing at and below the resonance. 
Therefore, a significant effect is predicted for disappearance at lower energies. 

Just such a lower energy z/ M disappearance experiment has been proposed [20] at the 
Spallation Neutron Source (SNS) being built at Oak Ridge. The neutrino source would 
be stopped 7r + 's, which undergo two-body decay to produce a monochromatic beam at 
30 MeV. In addition to the SNS source for stopped pions, there is the possibilty of a high- 
intensity "proton driver" at Fermilab which would also include stopped pions on its physics 
agenda. Detector distances at either site would be under 100 m from the pion source. Thus, 
the D/E is sufficiently small that only the LSND Sm 2 can effect neutrino flavor change. 

The amplitude for z/^-survival is given by eq. ([35)1. and the term oscillating with distance 
is sm 2 (5H D/2), with 5H given in eq. (|27jl. The effects predicted for the stopped-pion 
source at the SNS are shown in Fig. (0). The depletion of the beam due to substantial 
low-energy sterile-active mixing is considerable. For E <C E res , or for 8 near maximal, the 
oscillation length is insensitive to E ves . This explains the nearly common distance for the 
various minima in the figure. 

We note that the large z/^-depletion in Fig. (jUJ) is specific to the parameters we have 
chosen, and so should be interpreted as illustrative only. Smaller mixing leads to smaller 
depletion. If a v T component were added to the v a state, the //^-depletion may be less. 
Nevertheless, observable depletion of z/^'s from stopped pions is one of the more robust 
predictions of the brane-bulk model. 
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Figure 6: The muon-neutrino survival probability versus distance, for a monochromatic 
neutrino-beam energy of 30 MeV from stopped pions. The solid curves are parametrised 
by resonance energies 33, 100, 200 (identical for 300 and 400 MeV), in order of decreasing 
depletion. The dashed curve is the result with no bulk shortcut. The low-energy parameters 
for the 33 MeV resonance are chosen as in Fig. 3: sin 2 ^* = 0.01; sin 2 26 = 0.9; 5m 2 = 
0.7 eV 2 ; while the parameters for the other curves are chosen as in Fig. 4: sin 2 ^ = 0.1; 
sin 2 29 = 0.45; 5m 2 = 0.8 eV 2 . The vertical lines indicate the range of possible source to 
near/far detector distances. 



5 Further implications for astrophysics and cosmology 
5.1 BBN 



Successful big-bang nucleosynthesis puts severe constraints on the equilibration between 
active neutrinos and even a single sterile state, in the early Universe. The impact of 
active-sterile neutrino mixing on nucleosynthesis is quite complex. It has been discussed 
extensively, most recently in [21]. A popular idea to resurrect LSND in view of the BBN 
constraints is to introduce a lepton asymmetry, which gives effective masses to the active 
neutrinos in medium in the early universe, thus reducing the effective active-sterile matter 
mixing angles. 2 

The bulk shortcut effect will further differentiate the sterile and active neutrinos. It 

2 Another way to reconcile BBN with the existence of a sterile neutrino is to postulate a late-time phase 
transition [22]. With a reheating temperature of O(MeV) or less, the weak interaction simply does not 
have enough time to fully populate the neutrino modes. 
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also may provide an alternative to a lepton asymmetry. Consider natural expectations 
for the evolution of the brane metric: A higher density in the early universe will lead to 
greater gravitational attraction, and so to more brane buckles; and a larger temperature 
will increase thermal fluctuations of the brane. At the epoch of BBN, the temperature is 
ten orders of magnitude larger than today, and densities are thirty orders of magnitude 
larger than today. In the alternative metric (|18|) mentioned earlier, a higher density of 
scattering sites will increase the sterile neutrino scattering off from our brane. All these 
effects will increase the bulk shortcut parameter e and thus reduce the resonance energy 
(E res oc 1/y/e). If E ves is reduced to a temperature near enough the BBN temperature of 
~ 3 MeV, earlier oscillations will be suppressed. Our arguments that shortcuts had a larger 
e and therefore smaller E res in the earlier Universe are consistent with Ishihara's statement, 
that the magnitude of apparent causality violation on the brane increases with increasing 
matter density [7j. 

5.2 Sterile neutrino mass and dark matter 

The present scenario offers an effective mechanism for sterile neutrino production in the early 
Universe: when the temperature of the early universe drops below the resonant energy, the 
active-sterile mixing becomes maximal and active neutrinos are resonantly converted into 
sterile neutrinos. Afterwords, active neutrinos are re-populated via reactions maintaining 
thermodynamic equilibrium, until the active neutrinos decouple at energies around 1 MeV. 
The net effect, then, is to populate all neutrino modes, sterile and active. 

The occupation of states for a sterile neutrino with mass in the eV range impacts the 
effective total neutrino mass and number, which in turn impacts the connection ( "transfer 
function") between the cosmic microwave background anisotropy and today's large-scale 
structure. A similar impact of neutrino mass/energy obtains for measurements of galaxy 
bias stemming from galaxy-galaxy lensing, and for the large-scale power spectrum inferred 
from Lyman-alpha forest observations in the Sloan Digital Sky Survey. Specific conse- 
quences for cosmological evolution requires a detailed analysis (for recent works see |23j). 
Here we just mention that, as described in the previous subsection 15.11 higher temperatures 
and densities in the early universe may affect brane dynamics in a way to increase the 
effective shortcut parameter and to reduce the resonant energy. This effect may keep the 
sterile neutrino decoupled until the neutrino populations freeze out. 

Since the effective Sm1 SND ~ e V 2 is equal to ~ ^lsnd sm 2# on resonance, the true 
^ m LSND can De larger by 1/ sin 29. In the case of a small mixing-angle, the gain can be 
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very large. With a sufficiently small 8, the LSND sterile neutrino could play an important 
role as warm dark matter (WDM) with a mass of order keV. WDM has been proposed 
to solve the cuspy core problem of cold dark matter scenarios [21]. Sterile WDM has also 
been proposed to induce the observed high- velocities of radio pulsars j2H]. Unfortunately, 
it seems difficult to fit the actual LSND energy spectrum with a small mixing-angle 8, as 
the required small mixing angles inducing sharp peaks at the resonance energy, see Fig. El 

5.3 Supernova neutrinos 

Oscillations into sterile bulk neutrinos have interesting consequences for supernova neutri- 
nos. For example, supernova cooling may be accelerated due to the emission of Kaluza-Klein 
excitations of the sterile neutrino, resulting in constraints on a product involving the sterile- 
active neutrino mixing and the radius of the extra dimension and/or a delayed explosion 
process [26J. Since the efficiency of the bulk shortcut mechanism depends on the shortcut 
parameter e ~ (if) w ^ n only A being bounded from the radius of the extra dimension, 
such constraints can always be avoided by choosing a smaller A and a larger k. 

Furthermore, oscillations into sterile neutrinos would affect r-process nucleosynthesis, 
i.e. the rapid capture of neutrons on iron-sized seed nuclei, which is the prime candidate 
for the synthesis of nuclei heavier than iron. This process, which is believed to occur in 
type II supernovae, is suppressed by ^-capture on neutrons, which transforms the target 
neutrons into protons, forming stable a particles. It has been shown that ^-capture can 
be suppressed sufficiently if z/ e 's oscillate strongly into sterile neutrinos [2*7] . 

The bulk shortcut scenario will change this picture of r-process nucleosynthesis slightly. 
First, resonances, now involving matter effects and bulk effects, may find their energies 
shifted toward smaller values by the bulk shortcut. Moreover, supernova neutrinos with 
energies above the brane-bulk resonance will not experience any level-crossing when prop- 
agating out of the supernova, resulting in a cutoff of the active-sterile neutrino oscillation 
probability above E Tes . 

5.4 Horizon problem 

Finally, the sterile neutrino could couple more strongly to brane fields than the graviton, 
especially in the resonance region around E res . Thus, a solution to the horizon and homo- 
geneity problems, proposed in [SI El but based on gravitons, might turn out to be more 
effective if based on sterile neutrinos. Moreover, while bounds on the size of extra di- 
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mensions from precision measurements of the gravitational force law j2H] impose stringent 
constraints on the gravitational horizon, these bounds may not be valid for the extra dimen- 
sions felt by sterile neutrinos. The sterile neutrino horizon thus may be even larger than 
the gravitational horizon. Finally we stress that time dilation effects between the brane 
exit and the brane re-entry points can lead to causality violations which increase the sterile 
neutrino horizon. Such effects will be discussed in a forthcoming paper [29J. 




6 Discussion and Conclusions 

We have discussed active-sterile neutrino oscillation in an extra-dimensional brane world 
scenario. In such scenarios, sterile neutrinos paths may take shortcuts in the bulk, which 
imparts an energy dependence to the oscillation amplitude. Resonant enhancement of 
active-sterile neutrino mixing arises, parameterized by a shortcut parameter e = 5t/t. If the 
resonant energy lies in the range 30 MeV to 400 MeV, suitably chosen between the BUGEY 
and CDHS energies, then all neutrino oscillation data can be accommodated in a consistent 
3+1 neutrino framework. Such an energy range corresponds to e in the range 10" 18 — 10 -16 , 
and to brane fluctuations with a height to width ratio of ~ 10 -8 . The resonant energy might 
be identifiable in either the LSND spectral data and the muon neutrino disappearance from 
a stopped-pion source, or in the soon-to-appear MiniBooNE data. 

There are further interesting consequences for neutrino physics. We have mentioned 
that even the 2+2 model of active-sterile mixing can be resurrected with the brane-bulk 
resonance. Finally, we have sketched only briefly several interesting consequences for astro- 
physics and cosmology, consequences which remain to be worked out in detail. 
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